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ABSTRACT 

We carry out a numerical study of the effects of supernova-driven sliocks on galaxy 
formation at z = 9. These "cosmic explosions" can have a significant impact on galaxies 
forming nearby. We study such interactions in two key cases. In the first case in which 
the forming galaxy has already virialized, the impinging shock has only a small effect 
( ^ 1% of the gas is removed) and star formation continues relatively unimpeded. 
However, in the second case in which the nearby forming galaxy is at the more diffuse 
turn-around stage, a large fraction (« 70%) of the gas is stripped away from the host 
dark-matter halo and ejected into the intergalactic medium. As the time spent near 
turn-around is much longer than the interval from virialization to galaxy formation 
due to strong radiative losses, we expect the second case to be more representative of 
the majority of outflow-galaxy interactions. Thus SN-driven pregalactic outflows may 
be an efhcient mechanism for inhibiting the formation of neighbouring galaxies at high 
redshift. We briefly outline the possible cosmological consequences of this effect. 
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1 INTRODUCTION 

In currently favoured cosmological scenarios, the formation 
of structure is a hierarchical process, in which small initial 
subunits merge and accrete diffuse material to form larger 
structures of ever increasing mass. For instance, in the pop- 
ular Cold Dark Matter (hereafter CDM) cosmogony^, large 
numbers of objects with masses M = Kf'h'^ Mq collapse at 
~ 9 from 2-cr density fluctuations. In this collapse, the gas 
first falls inwards along with the dark matter, then shock- 
heats to the virial temperature (Tvir — 10** K; see, how- 
ever, Keres et al. 2004 who find that for low-mass objects 
a cold accretion mode can take place), condenses rapidly 
due to line cooling, and becomes self-gravitating. Subse- 
quently, massive stars form, synthesize heavy elements, and 
explode as supernovae (SNe) after ~ lO'^ yr, enriching the 
ambient medium. These processes may have been enhanced 
by primordial zero-metallicity stars (Population III stars) 



that may have been massive enough to i 


gnite powerful pair- 


production SNe (e.g.jHeEer & Wooslev 


f2005)), or generate 


laree 'seed' black holes (.Madau & Rocs 


2001). 



^ Throughout this paper we will assume a flat universe with total 
matter, vacuum, and baryonic densities in units of the critical 



density of Virn = 0.3, Ha 
a Hubble constant of Hq - 



-- 0.7, and Qi, = 0.05 respectively, and 
100 h km s-l Mpc-l, with h = 0.65. 



Around a redshift of 10-15, early sub-galactic stellar 
systems, perhaps aided by a population of accreting black 
holes in their nuclei, generated the ultraviolet radiation and 
mechanical energy that reheated and reionized most of the 
hydrogen in the universe. The latest analysis of the data 
from the WMAP ( Wilkinson Microwave Amsotropy Probe) 
satellite seems to suggest that thi s reionization took place at 

^ion — 

20l^° jKoeut et al.ll2003h . requiring the production 
of ionizing photons to have been extremely efficient in high- 
redshift galaxies (e.g. ICiardi. Ferrara fc Whit3 \2QQ^ ). 

The history of this crucial stage in cosmic structure for- 
mation depends both on the power spectrum of small-scale 
density fluctuations and on a complex network of poorly 
understood "feedback" processes. A generic consequence 
of all hierarchical scenarios is a large deposition of energy 
by SNe in the shallow potential wells of early sub-galactic 
systems. This may have two main effects, depending 
on the efficiency with which halo baryons can cool and 
fragment into clouds and then into massive stars: (a) the 
disruption of the host (dwarf) galaxy - the most violent 
version of stellar feedback - jLarsonI ITqtI iDekel fc Siu3 
119861: iMac Low fc Ferraral Il999l) : (b) the ejection of 
metal-enriched gas from the host galeixy and the subse- 
quent early pollution of t he int er galactic medium (IGM) 
jTeemark. Silk fc EvrarJ Il993l: [Cen fc OstrTkeJ Il999l: 
IScEmnEroiec^g^^a^u^l 1200 ll : l^bdEaTT^^^rafc Reed 
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[ 200 ll: IScannapieco. Ferrara fc Madaul 12002 

ISpringel fc Hernauisd |2003^ . Furthermore, the weU- 
estabhshed presence of heavy elements such as carbon, 
nitrogen and sihcon in the Lya forest clouds at z = 3 — 3.5 
provides strong evidence for such an e arly episode of 
pregalactic star form ation and outflows ijSongailal l200lt 
iPettini fc Pageill2004) . 

Similarly, stellar feedback and pregalactic outflows are 
likely to have efflciently inhibited high-redshift galaxy (and 
consequently star) formation. Because of the short cooling 
time of the gas in high-z small halos, strong feedback has 
been advocated in hierarchical clustering scenarios to pre- 
vent a "cooling catastrophe" in which too many baryons 
are converted into stars at very high redshifts. The required 
reduction of the stellar birthrate in halos with low escape 
velocities may naturally result from the heating and expul- 
sion of material by winds from massive stars and quasars, 
and multiple SN explosions. 

It is also well known that the radiative and mechan- 
ical energy deposited by massive stars and accreting black 
holes into the interstellar medium of protogalaxies may have 
a more global negative feedback on galaxy formation. The 
photoionizing background responsible for the reionization 
of the IGM will have two effects: a gas pressure increase, 
preventing accretion into low-mass halos, and the reduction 



los themselves ("Efstathiou 


'1992'; 'Thoul & 


Weinberg" 199^; 


[Navarro & Steinmetz 


1997 




Benson et al. 


2001 


n 


2Q02a b, 


l2003r. Furthermore, blast-v 


ves produced by mi 


i-quasars 



and protogalaxies propagate into the intergalactic space and 
may drive large portions of the IGM to a much higher 
adiabat than exp ected from photoionization (Voit 19J 



Madau. Ferrara fc Rees 2001: Theun s, Mo fc Schaye 2£ 
Cen fc Brvanll200lll . so as to quench the collapse of further 



galactic systems by raising the cosmological Jeans mass. 

In a complementary study, iBenson fc Madeiul ll2003ll 
have recently performed a detailed calculation of the effect 
of the increased gas pressure after this global early energy 
input into the IGM at the end of the cosmic dark ages. They 
have shown that preheating alone is not able to explain the 
sharp cut-off of the luminosity function at bright magnitudes 
and reproduce the observed abundance of satellite galaxies 
in the Local Group. It is worth noting that preheating con- 
sists of a homogeneous heat deposition in the IGM. There- 
fore, it affects galaxy formation in a glob al man ner, through 
an increase in the flltering mass |Gnedinll200(i l . 

In this paper we will be concerned with events of mul- 
tiple supernovae in protogalaxies and the mechanical in- 
teraction of their shocks with neighbouring overdense re- 
gions. Considerable recent theoretical effort has gone into 
understanding the long-range feedback due to UV/ionizing 
radiation. Relatively less attention has been paid to the 
problem of how galaxy formation in these objects is hin- 
dered by shocks from nearby protogalaxies. However, lim- 
ited analytic and numerical studies on a global cosmolog- 
ical scale have indicated that such interactions may have 
had an enormous im pact on the galaxy luminosity func- 
tion and its evolution ([S cannapieco^ Fer rara fc Madarl2002t 
iThacker. Scannaoieco .^^Davia2002^ !'^hile shock-cloud in- 
teractions within t he interstellar medium have b een simu- 
lated in detail Ce.g. lKlein. McKee fc Coellal lll994) 1. such re- 



lated cosmological interactions have been largely neglected. 
Thus it is important to study this process in detail. 

The structure of this paper is as follows: in i|21we outline 
the model used and the assumptions adopted in building up 
our model of "cosmic" stellar feedback; in i|3we describe the 
simulations carried out to examine the effects of SN explo- 
sions on a nearby halo, and in ij^we present our results and 
give a quantitative analysis. A summary is given in S|S| 



2 MODEL ASSUMPTIONS 

We want to study the effect of cosmic explosions produced by 
SNe in a given (source) halo impinging onto a neighbouring 
(target) object on its way to collapse into a luminous galaxy. 
To address this problem in a numerical context, we make 
use of the hydrodynamic code CLAWPACK'^, an explicit Eu- 
lerian, shock- capturing and upwind second-order Godunov 
(with TVD stability) method. 

Our investigation is centered on the impact on an im- 
pinging shock during two main stages of collapse of the tar- 
get halo: turn-around and virialization, both of which are 
described in further detail below. In both cases, we set the 
redshift of interaction to be z = 9, as a typical value during 
reionization, and the halo mass to M — Kfh'^ Mq, which 
is chosen both because of the ubiquity of such objects at 
z = Q (e.g. Scannapieco, Ferrara, fc Madau 2002) and be- 
cause they are able to quickly cool by atomic lines, regardless 
of the evolution of prim ordial H2 iHaiman. Rees. fc Loebl 
ll997tlCiardi et alJl200oll . 



2.1 Gas distribution in the target halo 

In constructing a model for the structural properties of the 
target halo, we neglect self-gravity of the baryons, which 
is motivated by the fact that ^i,/Q.m -C 1. Therefore we 
consider the gas in hydrostatic equilibrium inside the grav- 
itational potential of the dark matter halo. 



2.1.1 The virialized case 

To determine how the baryons are arranged in the dark mat- 
ter gravitational field, we must first determine the dark mat- 
ter distribution inside pre-galactic systems. Here we assume 
that dark matter halos at virialization equilibrium follow the 
universal (spherically averaged) density profile determined 
by Navarro, Frenk fc White (1997; hereafter NFW): 



p{r) 



cx (1 + cxY 



(1) 



where 



-Rvir 

pc is the critical density of the universe, iivir the (physical) 
virial radius of the system (at which the mean enclosed den- 
sity is 200 times the mean cosmic value = pc0.m{^ + z)^), 
c is the halo concentration parameter, 5^ = (200/3)c^/-F(c) 
the characteristic overdensity, and 

^ Conservation LAWs PACKage. Details can be found at 
http : //www. amath. Washington. edu/~claw/ 
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F{c) = ln(l + c) 



1 + c 



(2) 



The total mass of the halo within the virial radius and the 
mass within a radius r < R^ir are, respectively; 



and 



(1+2 



1'^ -Rvir 



(3) 



MM^f p(r')47rr"dr' = M:^^!^ . (4) 
Jo ^(c) 

Equation Q allows us to calculate a circular velocity 
squared 

where Vc = {GM / R^ivY^'^ is the value at the virial radius. 
The work done by the gravitational force on the escaping 
gas is equal to the variation of its kinetic energy, therefore 
the gas at the radius r is able to escape provided it has a 
velocity (squared) larger than: 



= (0 = 2 



GMjr') , _ 2 F[cx)+cxl{\ + cx) 



in which the NFW profile is extrapolated to infinity. The 
escape speed is maximum at the center of the halo, i'e(O) = 
2V'^clF{c). 

To proceed, we follow the algorithm described in the ap- 
pendix of NFW, and compute the concentration parameter 
(or, equally, the characteristic density contrast &c) of dark 
matter halos as a function of their mass for our adopted 
cosmological model. The algorithm assigns to each halo of 
mass M identified at redshift z a collapse redshift Zcoii, de- 
fined as the redshift at which half of the mass of the halo 
was first contained in pr ogenitors more mass ive than some 
fraction of the final mass ijLacev fc Colell993 ). The assump- 
tion that the characteristic density of a pre-galactic system 
is proportional to the critical density at the corresponding 
Zcoii implies: 



(7) 



Lower mass halos generally collapse at higher redshift, when 
the mean density of the universe is greater: for this reason, 
from the previous relation we deduce that low-mass systems, 
at any given time, are more centrally concentrated than mas- 
sive ones. As an example, we report in the following table the 
values of 2coii and c, obtained for different halos identified 
at z = 9. 



M[M0//i] 


2coll 


c 


10^ 


12.5 


4.9 


10» 


12.2 


4.8 


10» 


11.9 


4.7 



To study the impact of high-redshift cosmic explosions 
on objects at I + 2 ~ 10 in detail, we will assume a "typical" 
concentration parameter of c = 4.8. At these epochs, the 
dark matter halo of a subgalactic system is characterized by 
a virial radius: 



0.76 Mg/^ Q.-;^!'^ h'^ 



1 + 2 

10 



kpc. 



a circular velocity at virial radius: 

and a virial temperature: 
GM /imp 



(8) 



(9) 



where 
Ms = 



M 



108 Mq /l-l ' 

and fi is the mean molecular weight {fi = 0.59 for a fully ion- 
ized hydrogen/helium gas of primordial composition). The 
escape speed at the center is: 

l + 2\l/2 



.e(0) = 77A4''«J7^/«(i±£) 



km s ^. 



(11) 



Note, however, that high-resolution A'^-body simulations by 
iBullock et alJ l)200 J) seem to indicate that high-redshift ha- 
los are actually less concentrated than expected from the 
NFW prediction. If this is the case, we might slightly over- 
estimate the values of the escape speeds from pre-galactic 
systems. 

Having fixed the dark matter density distribution, we 
can establish the density profile of the baryons. If the gas 
collapses and virializes along with dark matter, it will be 
shock-heated to the viri al temperature and will settle down 
to an isothermal profile (|Makino. Sasaki fc Sutol|l998ll : 



Pb(r) = po e 
with 



-A [vm-vUr)]^ 



(12) 



A 



linip 



'^ik'J. vir 

Throughout this paper we make the assumption that the 
gas is distributed isothermally in virialized structures. The 
central density of the gas po is determined by the condition 
that the total baryonic mass fraction within the virial radius 
is equal to ilt/^m initially: 



Po 



(200/3) c3n6( 



Pc (1 + 2)3 JI^P{l+t)B/tdt 



44000 



(13) 



where B = 2c/F(c). At the virial radius pb{Rvii) = 
0.00144 po. 

2.1.2 The turn-around case 

Once a positive density perturbation has collapsed, it contin- 
ues to attract matter. Bound shells cease expanding with the 
Hubble flow, continue to turn-around, and consequently the 
mass of the perturbation grows by accretion of new material. 
This smooth accretion is generally called secondary infall in 
the literature, because particles fall onto an already formed 
object. In this case we make use of the results obtained by 
Bertschinger. (J,985) . Since gravity has no preferred scale, 
we can find self-similar solutions for secondary infall and 
accretion onto an overdense perturbation. The behaviour of 
this infall depends on the nature of the gas (collisional or 
collisionless) and on central boundary conditions, such as 
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the presence of a black hole in the core. In our situation, the 
coUisional (baryonic) matter moves in the potential well gen- 
erated by the dissipationless (dark- matter) component. The 
radius of the shell that is turning around increases with time 
as t**/^, while the mass within this radius increases as t ' . 
At large radii, both dark matter and baryonic components 
have the same density distribution, given by: 



pir) oc r ^^'^ 



(14) 



At smaller radii, the main difference between the two com- 
ponents is the existence of a spherical shock that forms in 
the baryonic gas due to the collision of the infalling gas with 
the self-generated core. Since this gas is dissipative, infalling 
fluid elements are decelerated by passage through this shock 
and by post-shock pressure gradients, and come to rest at 
the center. 

In this model, fluid motion is unchanged when param- 
eters are conveniently scaled as r ^ A = r/Rta,: 



Pb{r,t) = pigm/5(A) 
Mb(r,t) = ^TTpis^RtMiX) 



(15) 
(16) 

(17) 
(18) 



where v, pb, p, Mb are respectively the gas velocity, density, 
pressure and the mass at the distance r from the centre of 
the perturbation at the cosmological time t, while the same 
symbols with the "hat" represent the relative dimensionless 
quantities; i?ta is the turn-around radius and pigm is the 
mean baryonic density of the intergalactic medium at the 
redshift of interest (see equation l|2()|l'l. The shock is located 
at a fixed fraction of the turn-around radius, Aj, — 0.347 for a 
gas with 7 — 5/3, and propagates outwards according to the 
scaling rs oc t^^^ . The total mass of the over-dense region 
within the turn-around radius is calculated with Af (A — 
1) — 5.6 and with Qm instead of ^Ib in 11811 : 



3M 



4n 5.6 pci^n 



1/3 



(19) 



2.2 The intergalactic medium 

The intergalactic medium around the simulated halo is con- 
sidered to be homogeneous with a density pigm: 

Pigm = pc^b {1+zf = 1.88x10"^^ h^flb {1+zf g cm~^(20) 

corresponding to a numerical density of about 4 x 10"* cm~"^ 
at redshift z = 9. Furthermore, we assign to the IGM a 
temperature of T = 10* K, typical of a fully ionized medium. 
This value is explained by the observations of absorbing l inos 
of the IGM in spectra of distant quasars l)Rauchlll99?j l. It 
follows that the sound speed is Cigm ~ 15 km s~^. 



2.3 The shock wave 

To determine the properties of the impinging shock, we 
adopt a thin-shell approximation. In this model, we require 
an estimate of the shock Mach number M and the thickness 
of the cooled shell. If radiative losses are small within the 



hot-bubble driving the cosmic explosion, the shock velocity 
can be approximated by a Sedov- Taylor blast- wave solution: 

1/2 

d"^/', (21) 



V 37i-pig 

where E is the total energy of the shock, tfe = 0.2 (e.g. 
Scannapieco, Ferrara & Madau 2002) is the total-to-kinetic 
energy conversion efficiency, and K = 1.53 is a constant 
obtained from the exact solution for 7 — 5/3. As for the 
energy E, we assume that one supernova occurs f or every 
100 M0 of baryons that went into stars (see, e.g., ICibsonl 
1997), so the number A'^ of supernovae driving the bubble is: 

A = 5 X 10* Ms est (fi™ hy\ 

where e^f is the initial star formation efficiency. The total 
energy of the shock is then: 



E = 2xWekN erg. 



(22) 



where 2 x 10^^ erg is the energy of each supernova, 
taking into account also the contribution from stel- 
lar winds. Assuming an initial star formation efficiency 
egf = 0.1, a value that is broadly supported by ob- 
servations of high-2 I GM metals and th e high-z c osmic 



star formation rates l| Ci^d^^^ I2OOOI : iBarkanal |200_ . 
IScannapieco. Ferrara fc Madaul2002) . we can estimate N ~ 
5000 and thus E ~ 10^^ej;f2~^ erg for typical dwarf 
galaxy-sized objects of mass M ~ 10* M©. 

We estimate the distance d in equation (I21II . i.e. the 
radius of the shock wave, as the mean spacing between halos 
of a mass scale AI: 



d', 



( 3M 



\4-KpcO.-n 



j (1 + 2)"^ = 94M8/^(l + 2)-^ kpc. (23) 

From equations I20II . 1221 and 1231 1 we note that the value 
of Vs in equation H21|l is independent of redshift and halo 
mass. 

An alternative estimate relies on the fact that the mean 
number of halos within a distance r from another halo is: 



*(«; r) = 47rn / r'^ [1 + ^(r')] At' , 



(24) 



where n is the average number density of halos and ^ is 
their two-point correlation function. Both these quantities 
can be estimated from the classic analytical peaks formal- 
ism (Press & Schechter 1976, hereafter PS; Bond et al. 1991; 
Mo & White 1996). Note that the average number density in 
this equation is far less than (^d^)~^ as defined by equa- 
tion 1231 1. This is due to the fact that the mean separa- 
tion between halos is much larger than the average distance 
between a halo and its closest neighbor, because they are 
clustered together into groups of objects that only occupy a 
fraction of the cosmological volume. From equation 12411 . we 
can compute the probability for a nearest neighbour to be 
at a distance r as P(r) = 1 — exp[— 'I'(^ r)] and the mean 
distance to the nearest neighbour as: 



r— — (r) dr . 
dr 



(25) 



In this case, approximating the physical density n as 
(1 -I- z)~^ times the PS number density of -^^{M = 
10*Mq//i, z — 9), and making use of the matter power 



Dwarf Galaxy Suppression 5 



spectrum from Eisenstein & Hu (1998), with as = 0.87, we 
obtain a distance of 14 kpc. The similarity of this (power- 
spectrum and redshift dependent) estimate with our more 
straightforward approach, give us confidence that the dis- 
tance 9.4 kpc as given by equation 12311 is a reasonable esti- 
mate for interactions in which we are interested. 

We have to assess now the value of the second parame- 
ter of interest for the shock wave, the thickness of the shell 
Ad. Almost the totality of the gas swept up by the front is 
concentrated in a relatively thin layer and the density in- 
side the layer is constant and equal to the density behind 
the front surface. We have also to consider the fraction of 
baryons removed from the o utfiowing galaxy: according to 
iMori. Ferrara fc Madaul ll2002i) . we can take about the 60% 
of a system with the same mass scale M of interest. In con- 
clusion we obtain for the thickness Ad of the shell, from 
conservation of mass: 



density and the pressure pa behind the shock are 



4nig 



1.6 X 10 



cm ^ and ps 



Pigm Ws , respectively. 



Ad: 



^ Pigm _|_ 0.6 M flf 
3 Ps 47r Ps d2 



(26) 



where ps represents the post-shock density (i.e. the density 
inside the layer), achieved with the usual well-known "jump 
conditions" {ps — 4 pigm for supersonic shocks). 



3 SIMULATION RESULTS 

In both virialized and pre-virialized cases, we consider a 
computational grid made of 150 x 150 cells and set the 
boundary conditions as zero-order extrapolation, i.e. requir- 
ing that the flow of the matter runs only outward across all 
four borders. As mentioned above, in both cases we select a 
target halo mass of M = 10* ft" ^ Mo and redshift of z = 9 
as our fiducial model. Convergence tests with double spatial 
resolution of the grid (300 x 300 cells) have been made to 
ensure that the simulations provide reliable results: in all 
cases the answers are consistent. 



3.1 Virialized case 

Let us first consider the virialized case. For our fiducial pa- 
rameters the virial radius of the pre-galactic system, which is 
assumed to be spherically symmetric, is equal to 1.75 physi- 
cal kpc and its virial temperature is Tvir = 12, 490 K. In this 
case, the spatial resolution of the computational grid is 43.3 
pc. To determine the velocity at which the shock reaches the 
halo, we replace the distance d, in eq. (I23II with: 

d ^ d' = d - i?vir, (27) 

as we are interested in the initial point of the collision, when 
the front first reaches the overdense region. 
The simulation starts with the coordinates of the centre of 
the target halo placed at the point"^ {xo,yo) ~ (63,75), i.e. 
(2.7, 3.2) kpc, with a radius of ~ 40 cells. Both the pre- 
galactic system and the intergalactic medium are at rest, 
and the shock wave is modeled as a plane wave, at x = 20 
cells, moving from the left towards the right, with an initial 
velocity of about 214 km s~^. In this case, the numerical 

^ Hereafter we will indicate the coordinate of the cell along the 
horizontal axis as x and the coordinate along the vertical axis as 

y- 



and the temperature of the post-shock gas reaches a value 
of 6.2 X lO'^ K. 

Let t = he the time referred to by the initial condi- 
tions. The final time of the simulation is then t = 58.2 Myr, 
equal to about 3.6 times the shock crossing time of the halo, 
defined as: 



tc 



2-Rvir 



(28) 



Note that virialization is quite late in the evolution of the 
perturbation, and during our simulation time it is possible 
that star formation (and SNe) could occur inside this ob- 
ject. While, in this case, the dynamics of the system would 
become more complicated, we choose to ignore these efi'ects 
and study the virialized case as an idealization of the last 
possible time at which baryonic stripping could prevent star 
formation in the collapsing object. Indeed, due to the short 
time between virialization and star formation, it is much 
more likely that baryonic stripping will occur prior to viri- 
alization, as we discuss in more detail below. 

In Figure Q we show the density and temperature fields 
at three key stages. Initially the shock wave travels towards 
the centre, compressing the gas. When the shock first reaches 
the system, it encounters a density slightly greater than 30 
times that of the IGM. At this point, the velocity of the 
internal wave changes in direction and drops by a factor of 
approximately the square root of ratio of the halo pt and 
IGM Pigm densities (Bvc hkov fc Pikclneiiri975h : 



Vs,h 



(Pb/Pigm)^''^ 



(29) 



where Uj,, h and Vs^ igm are the front velocities inside the halo 
and in the IGM respectively, which gives ~ 38 km s~ . 

During the second stage, the shock velocity gradually 
decreases, according to equation 11211 as it runs into denser 
regions. In particular, the front reaches A4 ^ 1 close to the 
halo core and lapses into a sonic wave, mingling with the 
surrounding matter. At this time, the shock begins to sag 
and to form two symmetric extensions in which the baryons 
are accumulated. 

In the meantime the outfiow continues to propagate out- 
side the halo, with a velocity greater than the velocity inside 
the cloud. The intergalactic gas has now a post-shock pres- 
sure of ~ 1.4 X 10"^'' erg cm"'', more than 240 times the 
IGM pressure and about 6 times the baryonic pressure of 
the collapsed system at the virial radius. 

Finally, at a time t ~ 50 Myr the two flux sides join to- 
gether again behind the halo, "embracing" the entire system. 
After this time it is unlikely that the mass inside the virial 
radius will be ejected by momentum transfer. The dynamics 
of the interaction then becomes quite complicated, due to 
the presence of reflections of shock waves which scatter in 
the perturbation. 

In Figure |5| we show density and temperature pro- 
flies along the horizontal direction crossing the halo centre 
{y — 75) , at each of these three different stages (from top to 
bottom). These plots uncover the same trends Figure^ and 
indicate that after the interaction the gas structure of the 
halo is only weakly modified, especially in the core, where 
it is nearly unchanged. 
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3.2 Pre-virialized case 

Let us move now to the pre-virialized case, when the halo is 
at the turn-around stage. Here the radius of our fiducial sys- 
tem, with the same mass and at the same redshift, is given 
by equation 119L which yields Rta. ~ 5.7 kpc. Therefore, 
the Mach number of the shock wave given by equation I21II 
is larger than the previous value, and we obtain Ai = 37 
(i.e. a front expanding in the IGM with an initial velocity 
of roughly 565 km s^^ and with a crossing time of the halo 
tcross — 19.8 Myr) and Ad ~ 1.2 kpc. The post-shock tem- 
perature of the gas reaches the value of T ~ 4.3 x 10^ K. 
The new configuration of this problem consists of a spatial 
resolution of about 138 pc in the same 150 x 150 compu- 
tational grid. The centre of the overdense region is put at 
(49, 75) and its radius encloses 41 cells. The final time of the 
simulation is 133 Myr. 

The dynamics (Figure |^ of this interaction, by and 
large, follow the same sequence discussed in the virialized 
case. Nevertheless some important differences can be seen. 
Note that at this stage the halo is not yet isothermal and 



the solutions bv iBertschingerl il985l) expect values of pres- 
sure at the edge almost equal to those of the IGM, and a 
smoother decrease of the density with radius. 

It follows that at this epoch the shock will be able to 
more easily penetrate into the centre of the perturbation. 
Figure|3points out this fragility, showing that now the front 
is able to reach the halo core. For further details see also Fig- 
ure 2] representing the evolution in time (from top to bot- 
tom) of the density and temperature profiles along the hori- 
zontal cut crossing the halo centre {y — 75). The comparison 
between Figure and Figure |3 emphasizes a morphological 
difference of the system. 



4 QUANTITATIVE ANALYSIS 

As the primary aim of our investigation is to assess the im- 
pact of galactic outflows on neighbouring overdense regions, 
it is important to quantify the amount of matter that can 
be remov ed from the halo. We follow the two criteria pro- 
posed by IScannapieco. Ferrara fc Broadhursd ll2000l) to es- 




Figure 3. Snapshots of the dynamical evolution for the halo at the turn-around stage. Maps with logarithmic color scale are plotted for 
the density [g cm~^] (left) and temperature [K] (right) fields at three different times (t = 22.2, 44.4, 66.6 Myr). 
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Figure 4. Turn-around stage. Density (left) and temperature (riglit) profiles for a horizontal cut crossing the centre of the halo (y = 75), 
at times (from top to bottom) t = 22.2, 44.4, 66.6 Myr. Dashed lines show the profiles for the same cut at the initial time. 



timate this quantity. There are two mechanisms by which 
outflows from nearby objects can hinder the formation of a 
galaxy: in one case the gas is heated to a temperature greater 
than the virial temperature ("mechanical evaporation"), in 
the other case the shock transfers sufficient momentum to 
carry gas parcels out of the gravitational potential ("bary- 
onic stripping"). 

Figure |S] shows, for the virialized case, the total frac- 
tion of baryonic mass removed from the virial radius of the 
halo (filled circles), by each of these mechanisms. The tri- 
angles (diamonds) represent the evolution of the fraction of 
mass with T ^ Tvir, {v{r) ^ Ve{r)), integrated over the halo. 
The first condition describes the importance of the mechan- 
ical evaporation, whereas the second one refers to baryonic 
stripping. 

The physical interpretation of the behavior of the curves 
is as follows. The initial growth of the lines is due to the 
penetration of the shock into the outer layers of the halo, 
which are much less dense than the central core. Even at 
these large radii, however, the escape velocity is nearly 38 
km s~^ (see the lower curve), i.e. Ve ^ Vs 1291 1. Thus only a 



tiny fraction of the virialized gas is susceptible to baryonic 
stripping. 

As for the thermal contribution, while the gas on the 
outskirts of the perturbation is heated above Tvir, causing 
the initial rise in the upper curve, the short cooling times in 
the dense inner regions of the halo prevent T ^ Tvir heat- 
ing at later times. A heating plateau is soon achieved, after 
which additional baryonic heating becomes negligible. Thus 
at the end of the simulation {t = 58.2 Myr), only ~ 0.9% 
of the gas is affected by mechanical evaporation and only 
~ 0.7% is supplied with enough momentum to be stripped; 
hence, the impact of the cosmic blast is trifling. 

Figure|H]shows the analogous plots for the pre- virialized 
halo at the turn-around phase. In this case we find that 
galaxy formation is strongly affected, as a result of the 
stronger shock (the Mach number is larger than in the pre- 
vious case) and the more diffuse structure of the over-dense 
region, which make the system more vulnerable to the im- 
pinging blast. This is particularly evident in the case of bary- 
onic stripping, which increases remarkably (~ 69.9% at the 
end of the simulation, t = 133 Myr) in comparison with 
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Figure 5. Mass loss in tlie virialized case. Triangles represent the 
fraction of mass with T T^ii, while the diamonds points give 
the fraction of mass with v v^. Line with filled circles gives the 
total fraction of the baryonic mass removed from the virial radius 
of the halo by both mechanisms. 
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Figure 6. Mass loss in the halo at turn-around. Lines and points 
are as in Figure |21 Note the difference in scale on the y-axis. 

the virialized case. On the other hand the behaviour of the 
thermal contribution (curve with triangles) is only weakly 
enhanced with respect to the previous case, attaining a value 
of ~ 5.0% at the final output time. 



5 SUMMARY AND DISCUSSION 

In this paper we have studied the interaction of a cos- 
mic shock wave, generated by multiple supernovae occur- 
ring in a primordial galaxy, with a neighbouring halo ei- 



ther already virialized or at the turn-around stage. The 
main aim of the investigation is to assess if an intergalac- 
tic shock may be able to prevent the collapse of a fraction 
of dwarf galaxies at high redshifts, thus alleviating the so- 
cal led "cooling catastrophe" (fo r a review of feedback effects 
see lFerrara fc Salvaterral (|2004 )'). We have followed in detail 
the evolution of the target halo during and after the impact 
with the supersonic front by means of numerical simulations. 

The main results of this study can be summarized as 
follows. Suppression of dwarf galaxy formation can occur 
in the collapse stage at or before turn-around. In this case 
we find that about 70% of the gas can be stripped from 
the parent halo and ejected into the intergalactic medium. 
Most of the gas is unbound by the impinging momentum 
from the shock wave, which is sufficient to accelerate it 
to velocities larger than the local escape speed. The ef- 
fects of thermal evaporation of the gas (i.e. due to gas 
heating by the shock) are found to be negligible in com- 
parison to stripping, in agreement with t he conclusions of 
IScannapieco. Ferrara fc BroadhursO ll2000l) . 

The effect of the shock wave is much weaker if the halo 
is already in a virialized state because of the more compact 
configuration and consequent higher gas density, which en- 
hances radiative losses. The shock removes « 1% of the gas, 
and the forming galaxy is hardly affected by the interaction. 

As the time spent around turn-around is much longer 
than the interval from virialization to star formation due to 
strong radiative losses, we expect the pre-virialized case to 
be more representative of the majority of outflow-galaxy in- 
teractions. Thus SN-driven pregalactic outflows may be an 
efficient mechanism for inhibiting the formation of neigh- 
bouring galaxies at high redshift. 

Our numerical experiments have explored two typical 
situations among the possible shock - target galaxy interac- 
tions. However, they are by no means exhaustive. For exam- 
ple we have used mean values in order to choose the target 
galaxy (mass, formation redshift, density profile) and shock 
properties (energy, velocity at the interaction). Therefore 
a full exploration of the likely range for these parameters 
should be considered, at least for the more interesting case 
of pre-virialized halos, before a final conclusion on the rele- 
vance of these event can be drawn. Additional complications 
may arise from occurrences of multiple shocks impinging on 
the same halo, coming from different source galaxies. How- 
ever, given the relatively short duration of the interaction 
with the shock (100-200 Myr) the probability that more than 
one shock at a time is acting on the target galaxy seems un- 
likely, as this would require remarkable synchronization of 
the starbursts in different galaxies. 

The implications of the suppression effects might be 
particularly evident in the shape and evolution of the lu- 
minosity function. In particular, we expect a selective de- 
crease in the number of low-luminosity objects, due to their 
inability to collapse and form stars. Our results might also 
have implications for the number of satellites found around 
large galaxies, one of the major problems in current hierar- 
chical models of galaxy formation. The baryonic component 
of these satellites could be destroyed early on by interac- 
tions with shocks produced by the first cosmic episodes of 
star formation. Finally, stripping might drastically change 
the baryonic-to-dark matter ratios as a function of galactic 
mass. Although it is likely that SN explosions inside a given 
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galaxy will affect this observable, our study shows that im- 
pinging intergalactic shocks can play an equally important 
role. We plan to investigate these aspects in forthcoming 
work. 
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